o 

o 

(N 

oo 
(N 






cd 



X3 
O 

> 

in 

00 

o 
\o 
o 

-I— > 



I 

o 



X 



Theoretical study of interacting hole gas in p-doped bulk III-V semiconductors 
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We study the homogeneous interacting hole gas in p-doped bulk III-V semiconductors. The structure 
of the valence band is modelled by Luttinger's Hamiltonian in the spherical approximation, giving 
rise to heavy and light hole dispersion branches, and the Coulomb repulsion is taken into account 
via a self-consistent Hartree-Fock treatment. As a nontrivial feature of the model, the self-consistent 
solutions of the Hartree-Fock equations can be found in an almost purely analytical fashion, which 
is not the case for other types of effective spin-orbit coupling terms. In particular, the Coulomb 
interaction renormalizes the Fermi wave numbers for heavy and light holes. As a consequence, 
the ground state energy found in the self-consistent Hartree-Fock approach and the result from 
lowest-order perturbation theory do not agree. We discuss the consequences of our observations 
for ferromagnetic semiconductors, and for the possible observation of the spin-Hall effect in bulk 
p-doped semiconductors. Finally, we also investigate elementary properties of the dielectric function 
in such systems. 



I. INTRODUCTION 

Over the last years, effects of spin-orbit coupling 
in semiconductors have moved into the very focus of 
both experimental and theoretical solid-state research, 
mainly within the large and still rapidly growing field of 
spintronics^. In the p-type valence band of III-V zinc- 
blende semiconductors spin-orbit interaction is particu- 
larly strong. Important examples of p-doped semicon- 
ductor systems with itinerant charge carriers in the va- 
lence band include ferromagnetic semiconductors with 
Gai_a;Mna;As being the most intensively studied mate- 
rial, for an overview see Refs.^~^. In Gai_a;Mna;As and 
related systems, the substitutional Mn dopants form lo- 
cal moments with spin 5* = 5/2 from its five d-electrons, 
while they also act as acceptors providing holes in the 
valence band interacting with the local moments. This 
interaction between charge carriers and local spin mo- 
ments then leads, at low enough temperatures, to fer- 
romagnetic order, giving rise to the notion of carrier- 
induced ferromagnetism. So far. Curie temperatures as 
high as Tc ~ 160 . . . 170K have been observed^'^. There 
is a vast literature on the theoretical description of fer- 
romagnetic semiconductors taking into account realistic 
band structure parameterizations for the valence band, 
for an early key publication see Ref.^. However, what is 
most often neglected in the treatment of models for fer- 
romagnetic semiconductors is the Coulomb interaction 
among the holes. A (semi-)phenomenological way to ac- 
count for Coulomb repulsion is to introduce appropriate 
Fermi liquid parameters^. Exceptions to these heuristic 
approaches include numerical work based on dynamical 
mean field theory^, and a numerical Hartree-Fock study 
of a disordered two-band model, neglecting spin-orbit 
coupling^'^. For further dynamical-mean-field studies of 
models for ferromagnetic semiconductors not incorporat- 
ing Coulomb repulsion see Refs.^°~^^. 

Moreover, p-doped semiconductors have also attracted 



interest with respect to the recently predicted intrinsic 
spin Hall effect ^'*^^^; for a recent overview see also Ref.^^. 
In fact, the pioneering paper by Murakami, Nagaosa, and 
Zhang studies a p-doped bulk III-V semiconductor taking 
into account heavy and light hole bands around the F- 
point^"*. However, the Coulomb repulsion between holes 
was also neglected here. 

In summary, in the light of the above challenges and 
activities, it is certainly desirable to develop a deeper 
and possibly least partially analytical understanding of 
the effects of Coulomb interaction in p-doped semicon- 
ductors taking into account spin-orbit coupling. In the 
present work we study interacting holes in the valence 
band of a III-V semiconductor. The band structure is 
modelled by Luttinger's Hamiltonian in the spherical ap- 
proximation leading to heavy and light hole dispersion 
branches^*. The Coulomb repulsion between holes is 
treated via Hartree-Fock theory. As a nontrivial feature 
of the model, the self-conistent solutions of the Hartree- 
Fock equations can be found in an almost purely an- 
alytical fashion, which is not the case for other types 
of effective spin-orbit coupling terms. In particular, the 
Coulomb interaction renormalizes the Fermi wave num- 
bers for heavy and light holes. As a consequence, the 
ground state energy found in the self-consistent Hartree- 
Fock approach and the result from lowest-order pertur- 
bation theory do not agree. In other words, the self- 
consistent Hartree-Fock treatment contains contributions 
beyond lowest-order perturbation theory, which is a re- 
sult of the nontrivial band structure. We discuss the 
consequences of our observations for ferromagnetic semi- 
conductors, and for the possible observation of the spin- 
Hall effect in bulk p-doped semiconductors. Moreover, 
we also investigate elementary properties of the dielec- 
tric function in such systems. 

This paper is organized as follows. In section II 
we introduce the single-particle Hamiltonian and basic 
properties of the non-interacting system. In particu- 



lar, the structure of the single-particle eigenstates will 
be of importance for the Hartree-Fock study in section 
III. The self-consistent solution of the Hartree-Fock 
equations for Coulomb repulsion is presented in section 
III A. In section III B wc compare our findings for the 
three-dimensional hole gas with the situation in other 
generic semiconductor structures where spin-orbit cou- 
pling plays an important role. We then return to the 
three-dimensional hole gas and investigate its ground 
state energy and pair correlations functions in Hartree- 
Fock theory. In section HI E we also discuss elementary 
properties of the dielectric function within random phase 
approximation. We close with a discussion and outlook 
in section IV. 



II. THE NON-INTERACTING HOLE GAS 

A good approximative description of heavy and light 
hole states around the F-point in III-V zinc-blende semi- 
conductors is given by Luttingcr's Hamiltonian^* , 
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Here mo is the bare electron mass, p is the hole lat- 
tice momentum, and S are spin-3/2-operators, result- 
ing from adding the I — 1 orbital angular momentum 
to the s = 1/2 electron spin. The dimensionless Lut- 
tinger parameters 71 and 72 describe the valence band 
of the specific material with effects of spin-orbit coupling 
being included in 72. The above Hamiltonian is rotation- 
ally invariant and commutes with the helicity operator 
A = (fc • S)/k, where k = p/h is the hole wave vector. 

Thus, the eigenstates of (1) can be chosen to be eigen- 
states of the helicity operator, which grossly facilitates 
analytical calculations. The heavy holes correspond to 
A = ±3/2, while the light holes have A = ±1/2. For 
the dispersions of £ft/;(fc) of heavy and light hole states, 
respectively, one finds 
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where the masses to^ /; of heavy and light holes are given 

by 
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Well established values for the Luttinger parameters, 
among other band structure parameters, can be found in 
the literature^^. For example, for GaAs one has 71 w 7.0 
and 72 ~ 2.5 giving nih ~ 0.5m,o and mi w 0.08?Tio. 

The corresponding eigenstates of the Hamiltonian (1) 
are given by 



(r1fc,A) = -=|xA(fc)), 
V V 
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where V is the volume of the system. Using the conven- 
tional basis of eigenstates of S^- and introducing the usual 
parameterization k = k{cos(fsai'd^sanps\n'd,cosd) in 
terms of polar coordinates, the eigenspinors \x\{k)) of 
the helicity operator A = (fc • S)/k read explicitly 
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and the remaining eigenspinors |x-3/2(fc)), |X-i/2(fc)) 
can be obtained from the above ones by shifting -d i-^ 
TT — I?, ip i-^ (p + t:, corresponding to a spatial inversion 
k H^ —k. Note that \x±^{k)) is just a usual spin-coherent 
state of spin length S = 3/2 polarized along the direc- 
tion ±k/k. In what follows, we will also need the mutual 
overlaps squared between spinors which are given by 
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These expressions can be derived easily from Eqs. (5), (6) 
by putting one of the wave vectors along the z-direction 
and writing the resulting overlap squared in an explicitly 
rotationally invariant fashion as above. 

Let us now consider a non-interacting hole gas in an 
infinite system. Then the ground state is characterized 
by the Fermi wave numbers 



kh/i = 
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where £f is the Fermi energy. These wave numbers are 
related to the density n = N/V, N being the number of 
holes, via 
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(11) 



The kinetic energy per particle is straightforwardly ob- 
tained as 
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The above expression suggests to introduce an averaged 
mass TO by defining 

il+m^) (13) 

along with an averaged Fermi wave number 
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The kinetic energy per particle can then be rewritten as 
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which exactly resembles the familiar expression for the 
usual spin-1/2 electron gas^"'^^. In circumstances of 
Coulomb interaction between the holes, the above find- 
ing suggests to introduce a density parameter r^ and a 
Bohr radius as by defining 
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where e is the electron charge, and we have introduced a 
static dielectric constant Sr to account for screening from 
electrons in remote bands. Then the kinetic energy per 
particle can be rewritten as 
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where e'^/2aB£r is the Rydbcrg energy unit. Up to a 
slight difference in the prcfactor, the above expression is 
again completely analogous to the result for the usual 
electron gas^"'^^. However, as we shall see below, the 
exchange contribution from Coulomb interaction cannot 
be casted in a form immediately analogous to the spin- 
degenerate electron gas. 



to be translationally and rotationally invariant. Later on 
it will be specified to be the Coulomb repulsion. More- 
over, we assume a homogeneous neutralizing background 
ensuring charge neutrality and cancelling all direct (or 
Hartree) contributions from Hartree-Fock expressions. 

The eigenstates (4) of the single-particle Hamiltonian 
(1) solve the Hartree-Fock equations 



2 1,2 



erik){r^k,X) 



h-'k 

2toa 



(f|fc,A) 



Akr 



F(2^) 



-y 



(fk'{xy{k')\xx{k)) 



■V{\k^k'\)\xy{k')). 



(21) 



where T^(|fc|) is the Fourier transform of the interaction 
potential, and m\ stands for mu {mi) if A = ±3/2 (A — 
±1/2), a notation scheme which we will also use in the 
following. The Hartree-Fock eigenenergies are given by 
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To see that the eigenstates (4) solve the above Hartree- 
Fock equations, one can, again without loss of generality, 
take the wave vector k in Eq. (21) to point along the z- 
direction. Using the explicit parameterizations (5), (6) of 
the eigenspinors in terms of polar coordinates, one eas- 
ily sees that the integration over the azimuthal angle (p' 
ensures that the integral in Eq. (21) is indeed propor- 
tional to \xx{k)). This result holds for any interaction 
potential since |A: — fc'| is independent of ip' . The eigen- 
values (22) are then derived by performing the summa- 
tion over A' in Eq. (21) and using Eqs. (7)-(9). The 
first integral in Eq. (22) stems from the contributions 
with |A| = I A' I where as the second integral results from 
the cases |A| ^ |A'|. As explicitly shown in the nota- 
tion of Eq. (22), these eigenvalues arc functions of the 
integration boundaries qu^ qi arising in Eq. (21). In the 
presence of interactions, these quantities will in general 
not coincide with the Fermi wave numbers kh , ki of the 
non-interacting system, as we shall see below. 



III. THE INTERACTING HOLE GAS IN 

SELF-CONSISTENT HARTREE-FOCK 

APPROXIMATION 

We now consider an infinite system with a repulsive 
interaction between the holes which is naturally assumed 



A. Self-consistent Hartree-Fock solution for Coulomb 
interaction: renormalization of Fermi wave numbers 

Let us now specify the interaction to be the Coulomb 
repulsion, i.e. 
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The Hartree-Fock eigenenergies read explcitly 
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where the function h{x) is defined by 
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Note that h{l) is simply related to Riemann's (^-function, 
Ml) = (3/2)C(2)^7rV4. 

The two dispersions branches (24) coincide at zero 
wave vector, 
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for any values of qh, qi- However, when evaluated for the 
Fermi wave numbers kh, ki, they differ at the correspond- 
ing wave numbers, 
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Of course the Fermi energies for heavy and light holes 
have to be the same since otherwise a redistribution of oc- 
cupation numbers would take place. Therefore, in order 
to obtain a truly self-consistent solution to the Hartree- 
Fock equations, the Fermi wave numbers q^, qi have to 
be adjusted such that 



e/i iqh;qh,qi) ^ei iqi;qh,qi) 
under the constraint of a fixed density, 
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Thus, in fact just a single parameter, say qt, has to be de- 
termined numerically, which is technically a very simple 
task. Fig. 1 shows The ratios qh/i/khn of renormalized 
to unrenormalized Fermi wave numbers as a function of 
hole density for the III-V semiconductors GaAs, InAs, 



and InSb. The relevant parameters for these materials 
are summarized in table I. As seen from the figure, for 
realistic parameters one always has qu < kh and qi > ki, 
i.e. due to Coulomb interaction heavy hole states get 
depopulated in favor of light hole states. Moreover, the 
renormalization of Fermi wavcnumbcrs affects primarily 
the light hole wave number at low densities. The inset of 
Fig. 1 shows qh and qi as a function of density for GaAs. 
In Fig. 2 we have plotted the Hartree-Fock disper- 
sions E^,^{k;qh,qi) for GaAs at a hole density of n = 
5 • lO^^nm^^. The solid lines show the dispersion in- 
cluding Coulomb exchange for renormalized Fermi wave 
number q^/i , while the dashed lines represent the disper- 
sions of the free hole gas in the absence of interactions. 
As seen from Eqs. (24) the first derivative of the dis- 
persions £hn{k',qh,qi) with respect to k diverges both 
at k — qfi and k ~ qi giving rise to a vertical tan- 
gent at these points. In Fig. 2 these singularities are 
clearly pronounced for £h/i{k,qh,qi) at fc = q^/i while 
they are weaker and hardly visible in the plot at fc = qi/^. 
The fact that weak singularities occur in the derivative 
E^n {k;qh,qi) also at fc = qi/h is due to the mixing of 
heavy and light holes, i.e. the mutual overlap between 
heavy and light hole states at different wave vectors. 
Such an effect would be absent if one just had two spin- 
1/2-species of different mass, say electrons and muons, 
living in strictly different Hilbert spaces. 



B. Comparison with tiie two-dimensional electron 
gas and other systems 

As seen in Eq. (21), the eigenstates (4) of the non- 
interacting system provide solutions to the Hartree-Fock 
equations for a general pair interaction. This observation 
is familiar from the usual spin-1/2 electron gas with- 
out spin-orbit coupling^^'^i. If spin-orbit interaction is 
present, however, such a simple structure cannot be taken 
for granted, and the solutions to the Hartree-Fock equa- 
tions can in general become more complicated. As an ex- 
ample, consider a two-dimensional electron gas in a quan- 
tum well being subject to Rashba spin-orbit coupling^^. 
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where m is an effective band mass, and a is the Rashba 
parameter being tunable by an electric gate across the 
quantum well, a are the usual Pauli matrices describing 
the electron spin. We note that many-body effects in 
this type of system have recently attracted considerable 
interest2'^~29_ 'pj^g above Hamiltonian has two energy 
branches, 
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with eigenstates 
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where A is the area of the system. Now consider an arbi- 
trary pair interaction. Since the total Hamiltonian is still 
translationally invariant, the solutions to the Hartree- 
Fock equations can always be chosen to have a form 
similar as above, i.e. a plane wave factor times a two- 
component spinor. However, it is easy to see that the 
eigenstates (32) do not provide solutions to the Hartree- 
Fock equations for a general interaction potential. This 
observation is due to the fact that the angular integration 
in this two-dimensional case is different from the three- 
dimensional situation of Eq. (21). Only for a pure con- 
tact interaction (having a constant Fourier transform), 
the eigenstates (32) solve the Hartree-Fock equations. 

The same conclusions apply to bulk valence-band elec- 
trons being subject to the three-dimensional Dresselhaus 
spin-orbit coupling term^*^, and to heavy holes in asym- 
metric quantum wells'^^''^^. In both cases, the effective 
spin-orbit interaction is trilincar in the particle momen- 
tum. 

The situation of the two-dimensional electron gas be- 
comes even more complicated if also Dresselhaus spin 
orbit coupling is considered which reads in its two- 
dimensional approximation^"^ 
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with a coupling parameter /3. However, in the case when 
the Rashba parameter is of equal magnitude as the Dres- 
selhaus parameter, a = ±/3, the corresponding eigen- 
states of the non-interacting system solve the Hartree- 
Fock equations for an arbitrary interaction potential. 
This is due to the additional conserved quantity aris- 
ing at this point'^'* which cancels the effects of spin-orbit 
coupling in many respects^**'^^. 



C. Total ground state energy 

Let us come back to the case heavy and light holes 
interacting via Coulomb repulsion in the valence band of 
bulk HI-V semiconductors. With the renormalized Fermi 
wave numbers qh, qi, the total kinetic energy per particle 
in the ground state reads 
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with the density n given by Eq. (29), and for the total 
exchange energy per particle one finds 
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resulting in a total energy per hole 
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To obtain the corresponding results for the unrenormal- 
ized Fermi wave numbers, one just has to replace in 
the above expressions q^n with ki^n. Note that the un- 
renormalized expression Etot{kh,ki) is the equivalent to 
first-order perturbation theory in the Coulomb repulsion 
where one just computes the expectation value of the 
interaction with respect to the ground state of the non- 
interacting system characterized by the Fermi wave num- 
bers kh/i- 

Fig. 3 shows the ground state energy per parti- 
cle Etot{qh,qi)/N from the self-consistent Hartree-Fock 
treatment and the result Etot{kh,ki)/N from lowest- 
order perturbation theory as a function of the density for 
GaAs. As seen in the figure, it is always Etot{qh,qi) < 
Etot{kh,ki), i.e. the self-consistent Hartree-Fock ap- 
proach gives the lower ground state energy. This is clear 
since the ground state obtained from a self-consistent so- 
lution of the Hartree-Fock equations is by construction 
the Slater determinant of lowest energy in the Hilbert 
space of the many-particle system. Thus, any other 
Slater determinant state has to have a higher energy ex- 
pectation value. In fact, the renormalized Fermi wave 
numbers can alternatively be obtained by minimizing 
Etotiqh, qi)/N with respect to qh, qi under the constraint 
of a fixed density n. In the absence of interactions the 
minimization of Ekin{qh,qi)/N immediately reproduces 
the results of section II (as it has to be), whereas for 
Etot (qh 1 qi)/N one ends up with a coupled system of poly- 
nomial equations which does not seem to allow for an 
explicit analytical solution. 

The hole densities usually occurring in realistic sam- 
ples lie at n = 0.01 . . . l.Onm^^. At such densities, the 
difference between Etot{qh,qi)/N and Etot{kh,ki)/N is 
indeed very small (cf. Fig. 3), and the density parameter 
rg as defined in Eq. (18,(19) is of order unity, giving con- 
fidence to the validity of the Hartree-Fock approach^°'^^. 
At smaller densities like n < O.OOlnm""^, Etot{qh,qi)/N 
and Etot{kh,ki)/N differ appreciably. However, at these 
densities Hartree-Fock theory cannot be expected to give 
accurate results. On the other hand, it is common in 
many-body perturbation theory to refer to all contri- 
butions to the ground state energy beyond the lowest- 
order exchange term as correlation contributions^"'^^. 
In this sense the difference between Etot{qh,qi)/N and 
Etot{kh,ki)/N (resulting from the renormalization of 
Fermi momenta) can be viewed as a correlation effect. 

The inset of Fig. 3 shows the same data as the main 
panel, but as a function of the density parameter r^. The 
minimum of the Hartree-Fock ground state energy lies at 
about Ts ~ 5, similarly to the case of the usual spin- 1/2 
electron gas^'^'^"'^. Note that the maximum difference be- 
tween Etotiqh, qi)/N and Etot{kh,ki)/N is also achieved 
around this value. 



D. Pair distribution function 

It is instructive to also investigate the pair distribution 
function g(r) defined by 
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where I, J label the particles in the system and (•) de- 
notes the expectation value within the ground state. The 
ground state obtained from self-consistent Hartree-Fock 
theory is a single Slater determinant. Here the pair dis- 
tribution function can be formulated as 
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where gf^f^ igft) s-re the exchange contributions from 
heavy (light) hole states only, whereas gf^^ stems from 
exchange between heavy and light holes. It is straight- 
forward to calculate these contributions explicitly using 
Eqs. (7)-(9). The results can be formulated as 



9hl{r) 




^ihiqhr)f + ^ihiq,r)f 



+ fih{q,r)f 
"mi\^ I 9 



(39) 



-^'^(^)-(f)^U(^^(-^^^-S(^^(-» 



+ '^i{h{qir)f]. 



(40) 



9^nr)^i^f\^T 'iihiqnr))^.iq.r)) 
\ m / \ m / 



ifi J \ 16 

27 
-y(/2(g.r))(/2(g,r)) 

.^{I,{q,,r)){h{qir)) 



where we have defined 
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Note that 5(0) = 3/4, corresponding to a fcrmionic gas 
with four spin components. In the absence of spin-orbit 
coupling, qh — qi , the contributions involving I2 , I3 can- 
cel and one obtains the well-known exchange terms of 
the usual electron gas^^'^"'^. Fig. 4 shows the pair distri- 
bution function in GaAs for three different densities. At 
high enough hole density one can see Friedel-type modu- 
lations of g{r) whose period is essentially given by twice 
the heavy-hole Fermi wave number q^- 



2 a;2 



2 a;3 



E. The dielectric function in random phase 
approximation 

Within random phase approximation (RPA), the di- 
electric function is given by^°'^^ 



e«^-^(fc,c^) = l-y(fc)xo(fc,w). 
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where xoik, to) is the susceptibility of the non- interacting 
system. Its real part has the form 
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Here /(fc, Ai) are Fermi functions, and to obtain Eq. (46) 
from Eq. (47) we have used elementary properties of the 
spinor overlaps and the dispersion relations of the non- 
interacting system. In particular, in the static limit one 
has 
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By construction, xa{k,(-u) is entirely determined by the 
properties of the non-interacting system. In particular, 
at zero temperature, the integration boundaries in the 
above expressions are given by the unrcnormalized Fermi 
wave numbers kh, h. However, in order to be consistent 
with the self-consistent Hartree-Fock treatment, one may 
use the renormalized wave numbers qh, qi instead. As 
seen above, at high enough densities, the difference is 
negligible. 

For early work on dielectric response in zero-gap semi- 
conductors we refer to Refs^®'^^. An evaluation of the 
static expression (48) for the case of the two-dimensional 
electron gas with Rashba spin-orbit interaction has been 
given by Chen and Raikh^'^. Their findings have recently 
been challenged by Pletyukhov and Gritsev^*. The main 
technical obstacle there is posed by non-elementary inte- 
grals. In the present case of the three-dimensional hole 
gas, however, the occuring integrations are mostly ele- 
mentary but often very tedious. 

Analogously as for the three-dimensional elelctron gas 
the static dielectric function for the hole system at zero 
temperature can be formulated as 



.RPA 



(fc,0) = l + 






L 



k 



k 
2^1 



where 
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is the usual Thomas- Fermi screening wave number. The 
function L is the analogue of the well-known Lindhard 
correction for the electron gas. The explicit form of L 
for the hole system, however, is rather lengthy and te- 
dious and shall not be given here. For long wavelength, 
fc — ^ 0, L approaches unity. In this limit terms mixing 
heavy and light holes do not contribute to the dielectric 
function. A similar finding is valid for the plasma fre- 
quency which characterizes collective excitations at zero 
wave vector^'^'^"'^. Here contributions mixing heavy and 
light holes are also absent, and the plasma frequency ujp 
is given by 
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(51) 



where Uh, ni are the densities of heavy and light holes, 
respectively. 



IV. DISCUSSION AND OUTLOOK 

We have studied the homogeneous interacting hole 
gas in p-doped bulk III-V semiconductors modelled by 
Luttinger's Hamiltonian in the spherical approximation. 
The Coulomb repulsion is taken into account via a self- 
consistent Hartrcc-Fock treatment. As a nontrivial fea- 
ture of the model, the self-consistent solutions of the 
Hartree-Fock equations can be found in an almost purely 
analytical fashion. As we have discussed in detail in sec- 
tion IIIB this is not the case for other types of effec- 
tive spin-orbit coupling terms. As an important quali- 
tative feature, the Coulomb repulsion renormalizes the 
Fermi wave numbers for heavy and light holes: The in- 
teraction leads to a redistribution of occupation numbers 
from heavy holes to light holes compared to the non- 
interacting case. As a consequence, the ground state en- 
ergy found in the self-consistent Hartree-Fock approach 
and the result from lowest-order perturbation theory dif- 
fer from each other. By construction, the self-consistent 
Hartree-Fock result gives the lower ground state energy. 

The three-dimensional III-V semiconductor hole gas 
in particularly relevant for ferromagnetic semiconductors 
which are usually p-doped materials. In the theoretical 
description of these materials, the interaction between 
itinerant charge carriers is most often neglected^"'*, or ab- 
sorbed in effective Fermi liquid parameters*. The typical 
hole densities in Mn-doped GaAs, the most prominent 
and best-studied ferromagnetic semiconductor, arc of or- 
der n « O.lnm^'^. For such carrier concentrations, the 



density parameter rg is of order unity, giving confidence 
to the validity of the Hartree-Fock approach^"'^^. More- 
over, as seen in the present investigation, the renormal- 
ization of Fermi wave numbers is negligible at such den- 
sities, i.e. the interacting ground state in Hartree-Fock 
approximation and the non-interacting ground state are 
practically the same. In this sense, the abovementioned 
models for ferromagnetic semiconductors neglecting the 
Coulomb interaction are supported by the present study. 
However, the single-particle Hamiltonian used there is 
a simplified one which does not take into account the 
split-off band. A more complete description of the va- 
lence band including these states is given by the six-band 
Kohn-Luttinger modeP* as used in Refs.^'*. In fact, the 
influence of the split-off band is known to be important 
for the stability of the ferromagnetic order"^^. However, 
for the full six-band Kohn-Luttinger Hamiltonian mainly 
analytical progress like in the present work is certainly 
not possible, and one would need to resort to more com- 
plicated numerics 

During the last years various predictions and subse- 
quent experiments regarding spin-Hall transport in semi- 
conductor systems have attracted a very remarkable deal 
of interest^^"^^; for a recent overview see also Ref.^^. The 
first work opening the field of intrinsic spin-Hall effect 
was a paper by Murakami, Nagaosa, and Zhang who 
considered a p-type bulk III-V semiconductor^^. The 
single-particle Hamiltonian used in Rcf.^^ is the same as 
here with the Coulomb repulsion between the holes be- 
ing neglected^'*. For a disorder-free system, the spin-Hall 
conductivity is given by^^'"**^ 
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where the direction of the spin current, its polarization 
direction, and the direction of the electric field are mutu- 
ally perpendicular. In the above expression we have used 
the renormalized Fermi wave numbers qh, qi- Figure 5 
shows a spin-Hall conductivity as a function of hole den- 
sity for GaAs both for renormalized and unrcnormalized 
Fermi wave numbers. To facilitate the comparison to the 
usual charge conductivity, we have converted the spin- 
Hall conductivity to units of charge transport by multi- 
plying with a factor of e/ti. As shown in the figure, at 
densities n ^ O.Olnm^'' typical for realistic samples, the 
difference between the case of renormalized and unrcnor- 
malized wave numbers is negligible. Appreciable discrep- 
ancies occur only at small densities, where the validity of 
the Hartree-Fock treatment becomes questionable any- 
way. 

Finally, we hope that the present study will initiate 
further investigations on interacting III-V semiconductor 
hole systems. One possible direction is to perform (pre- 
sumably numerical) Hartree-Fock calculations for more 
complex band structure models as mentioned above. An- 
other obvious goal for future studies is to investigate 
many-body effects beyond the Hartree-Fock level. 
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FIG. 1. The ratios of renormalized and unrenormalized 
Fermi wave numbers, g^/j and fc^/j, respectively, as a func- 
tion of hole density for different III-V semiconductors. The 
renormalization of Fermi wave numbers affects primarily the 
light hole wave number at low densities. The inset shows qn 
and g; as a function of density for GaAs. 
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TABLE I. The Luttinger parameters 71, 72, effective hole 
masses m^/j, and static dielectric constants Sr for the III-V 
semiconductors GaAs, InAs, and InSb. 
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FIG. 2. The Hartree-Fock dispersions 
GaAs at a hole density of n = 5 • 10~^nm~'^. The solid lines 
show the dispersion including Coulomb exchange for renor- 
malized Fermi wave number Qh/i, while the dashed lines rep- 
resent the dispersions of the free hole gas in the absence of 
interactions. The singularities in £^,^{k;qh,qi) at A: = g^/j 
are clearly pronounced while the singularities at at A; = qi/h 
are weaker and hardly visible in the plot. 
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FIG. 3. The ground state energy per particle Etot{qh, qi)/N 

from the self-consistent Hartree-Fock treatment (with renor- 

malized Fermi wave numbers) and the result Etot(kh,ki)/N 

from lowest-order perturbation theory (with unrenormalized 

Fermi wave numbers) as a function of the density for GaAs. 

The inset shows the same data as the main panel but as a 

function of the density parameter Ts- 
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FIG. 4. The pair distribution function g{j-) for GaAs at 
three different densities. 
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FIG. 5. The spin-Hall conductivity a (converted to units 
of charge transport) as a function of hole density n for GaAs. 
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